Hohenberg-Kohn theorem for the lowest-energy resonance of unbound systems 
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We show that under well-defined conditions the Hohenberg-Kohn theorem (HKT) can be extended 
to the lowest-energy resonance of unbound systems. Using the Gel'fand Levitan theorem, the 
extended version of the HKT can also be applied to systems that support a finite number of bound 
states. The extended version of the HKT provides an adequate framework to carry out DFT 
calculations of negative electron affinities. 

PACS numbers: 



Most ground-state properties of electronic systems can 
now be calculated from first principles via Density Func- 
tional Theory (DFT) [1]. When the ground state of an TV- 
electron system is not bound, strict application of DFT 
with the exact exchange-correlation functionial should 
yield results that are identical to those of the (TV — M)- 
electron system, where TV — M is the maximum number 
of electrons that the external potential can bind. But in 
such cases, rather than the absolute ground state one is 
often interested in resonant states (long-lived metastable 
states), even if not eigenstates of the Hamiltonian. It is 
the lowest-energy resonance (LER) that plays the role of 
the ground state in the sense that, during its lifetime, it 
best represents the physical state of the TV-electron sys- 
tem. In fact, the LER is the localized "ground state" of 
the non-hermitian operator that is obtained by complex- 
scaling the coordinates of the original TV-electron Hamil- 
tonian by an appropriatephase factor (for complex scal- 
ing techniques, see refs.|2|, |3|, lj|). The energy £ LER and 
inverse lifetime r LER of the LER are given respectively 
by the real and imaginary parts of its complex eigen- 
value -E L er- One can associate a complex density ng(r) 
to it, and intuition suggests that the Hohenberg-Kohn 
theorem (HKT) that provides the foundation of DFT [5| 
can be extended along the following lines: the complex 
density ng(r) associated with the LER uniquely deter- 
mines the 0-scaled external potential, and all properties 
of the LER (in particular £ LER and r LER ) are therefore 
functionals of ng(r). History teaches us to be watch- 
ful, however, since the HKT has proven elusive when at- 
tempting to depart from the ground state. It has been 
shown to hold for the lowest-energy state of any given 
symmetry 0, 0] , with the resulting functionals depend- 
ing on the particular quantum numbers corresponding 
to each symmetry. More importantly, the lack of an 
HKT for excited states was recently demonstrated 0, i.e. 
excited-state densities, in general, do not uniquely deter- 
mine the external potential. We are faced here with a 
rather different problem, since the LER is an eigenstate 
of the complex-scaled Hamiltonian rather than the un- 
sealed one. The simplest example is the 3 P resonance of 



H~ (the lowest-energy state for that symmetry Q), for 
which ground-state DFT would predict its energy, had 
nature not made it unbound. 

In spite of the fact that the HKT holds only for the 
ground state [l(J, attempts have been carried out to use 
time-indepedent DFT for the calculation of excited en- 
ergy levels , and even very recently for resonances that 
can be regarded as excited states where the widths (in- 
verse lifetimes) are not equal to zero, but still small 
Whereas our aim is similar in spirit to that of ref. [l2| , we 
focus our attention on the LER of unbound systems. 

Complex-coordinate scaling is a well-developed tech- 
nique to characterize resonant states: upon multiplying 
all electron coordinates of the Hamiltonian by a phase 
factor e , the complex-scaled Hamiltonian Hg has right 
and left eigenvectors, denoted by \^$) and ($?g\ respec- 
tively, at: (a) bound states of the original (9 = 0, her- 
mitian) Hamiltonian, corresponding to the same energy 
eigenvalues; (b) continuum states of the original Hamil- 
tonian; the respective eigenvalues are rotated into the 
lower-half of the complex energy plane by an angle of 26; 
and (c) resonant states, with ^-independent eigenvalues. 
The complex-scaled resonance eigenfunctions are expo- 
nentially localized in the interaction region, whereas the 
continuum eigenfunctions almost vanish there [1 31 ] . We 
emphasize that the variational theorem for non-hermitian 
quantum mechanics has been developed only for reso- 
nances and not for the continua[l5j]. 

The original Hohenberg-Kohn proof Q is based on the 
minimum principle for the ground-state energy and on 
the fact that the TV-electron density operator n(r) — 
£i=i ^( r ~^i) couples linearly with the external potential 
^ext (r) , i.e. that one can always write a static TV-electron 
Hamiltonian as: 

H = f + V ee + J drh{r)v cxt (r) , (1) 

where T is the TV-electron kinetic-energy operator, and 
V ee is the electron-electron repulsion. This linear cou- 
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pling is of course maintained upon scaling (z = re l6 ), 

He =Tg + V ee ,e + / dzn(z)v cxt (z) , (2) 



where Tg 



-2i6 



T, V a 



-iO 



V ee for Coulomb in- 



teractions, and n(z) = exp(— z36>)n(r)[14]; but the mini- 
mum principle no longer holds: the complex variational 
principle [15| guarantees stationarity at all the resonant 
eigenf unctions of He , but not minimality at any of them. 
The main result of this paper is the realization that mini- 
mality at the LER is generally true for unbound systems, 
and, as a consequence, a practical analog of the HKT can 
be established. 

Minimality at the LER. Start with a simple case to 
motivate our statements. Consider a single electron mov- 
ing in the one-dimensional potential 



v C xt(x) = Qx 2 -a^je 0x2 ,x>0, 



(3) 



(vcxt — > oo for x < 0), and choose a and (3 so that i> ex t 
has no bound states. Fig.l shows this potential, along 
with the complex energies and magnitude squared of five 
resonance wavefunctions obtained via complex-scaling for 
an appropriate choice of a and (3. The LER energy for 
such choice is £ LER = Re(^f ER \H e \^ e LEU ) = 0.62a.u. We 
ask whether the expectation value of He for an arbi- 
trary trial square-integrable function $ tr iai can have a 
real part that is less than £ LER . Choose for example 
<J> trial (a;) = Cxe"> x2 and set C so that $f iiaJ = <J> trial (xe ie ) 
is properly normalized. We show in Fig. 2 the energy 
Atrial = R e ($ t *rli 1^0 1 ^taiai) as a function of 7 and note 
that it is above £ LER for all 7. According to the bounds 
derived by Davidson et. al. [Lsl ] for resonance positions 
and widths, there is no reason to expect this to be always 
the case, since the most one can say about the exact com- 
plex eigenvalue at a resonance is that it lies within a circle 
of radius determined by the complex variance associated 
with the trial wavefunction. But the LER of unbound 
systems is a special resonance since no 0-independent 
eigenvalues of He exist below it. A local minimum at 
the LER must also be a global one. Nothing guarantees, 
however, that the energy at the LER is a local minimum, 
rather than a local maximum, or saddle point. 

We now argue that the result observed in Fig. 2 for our 
test example is in fact usually the case, also for N elec- 
trons, and discuss the plausibility of this statement for a 
special subset of trial functions: those that are localized 
in the region where the resonance wavefunctions, and in 
particular the LER, have a high amplitude (see our com- 
ment above on the localization of the resonances in the 
interaction region, whereas the rotating continuum states 
are not localized). To be specific, restrict the discussion 
to trial functions that satisfy: 



\(KL\^)\ 2 >1 



(4) 



The same condition was employed before in ref.[l8| to 
derive upper and lower bounds for resonances. In such 
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FIG. 1: Top panel: potential of Eq.(|3]) (thick line) and mag- 
nitude squared of 5 resonance wavefunctions (a = 0.8 and 
/3 = 0.1). Bottom panel: Energy spectrum of He for 6 — 0.5. 
There are no bound states. The continuum branch cut is ro- 
tated by —28 with respect to the real axis, and 5 resonances 
are clearly exposed. 
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FIG. 2: The real part of the trial energy is above the LER 
energy for all values of 7 (see text). 



cases, Atrial may be expanded in terms of only resonance 
states. There are no bound states to populate, and the 
overlap with 0-scaled continuum eigenstates is negligible. 
(On the fact that the resonances serve as an almost com- 
plete set, see Ref . [l9| ) . Then, to a good approximation, 



E. 



trial J - / LER 



1+ J2 (n 8 \< M ) 2 (Ek/E^ R -i) 

fc^LER 



(5) 

where the ^ are resonance eigenfunctions of He with 
complex eigenvalues Ek — £*; — f Tfc. We conclude that if 



£ [Re(*f|< mI ) 2 (£ fe -£ LER ) 

fc^LER 

1 



then 



|i m (*f|$f rlal ) 2 (r fc -r LER ) >o ,(6) 



Atrial ^ ^LER. (7) 
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of the Hohenberg-Kohn theorem follows. Two potentials 
Ugxt and Vg^t that do not support any bound state, and 
differ by more than a constant, cannot yield the same 
LER-density n e (r) = (*£LK r e ie )l*L E R>- To see this, 
assume that the two potentials could in fact give rise to 
the same LER-density: 



<*Z&|TKre*)|*£k> = <¥ZSWre")|*SL) 



(9) 



where ^le R 



is the LER-eigenstate of the Hamiltonian 



of Eq.([2]) with external potential 



Using *y R as 



FIG. 3: Top pannel: Total overlap S of a trial function 



a trial function to estimate the lowest-energy eigen- 
value of we get by virtue of Eqs.© and © that 



with the resonance eigenfuctions Re^Lal^e + V eefi \^L) > R-e(**^|T e + V ee> |*ler) 



"[2] I 



,[2] 



of H 6 : S = ELi(*» I*') • The expansion is adequate in 
the shaded region. Deviations from S = 1 for large 7 diminish 
when broader resonances are included in the expansion [19(. 
Bottom panel: left-hand-side of Eq.© (thick line), along with 
the contribution to it from the real part of the overlaps (solid) 
and the imaginary part (dotted). 



The inequality of Eq.© would be obviously true if we 
were dealing with bound states, since then all the inverse 
lifetimes would be zero, all the real parts of the over- 
laps would be positive, and £k — f LBR > for all k, by 
definition of the LER. 

Eq.©, and therefore Eq.([7]), also hold true for res- 
onances because: (1) Normalization of the trial func- 
tion, Efc Re (*fc e | $ ?riai) 2 = !, together with the con- 
dition given by Eq.Q and the fact that — £ LER is 
positive for all k > 1 (fc = 1 denotes the LER), lead 



LER / 
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to E^iMn ^ 



y(£ k - £ LER ) > 0; and (2) Nor- 



malization also requires ^2 k Im(^^, £ '|$f r . al ) 2 = 0, so if 
T/t^i — F LER is a smooth function of k, as is usually the 
case, then the second term of Eq.© is expected to be 
smaller than the first one. 

We illustrate all this in Fig. 3 for our one-electron 
toy example. The two conditions discussed above are 
shown to hold in the region where an expansion in terms 
of resonance eigenstates is adequate (shaded region in 
the figure). Based on the result of Fig. 2 showing that 
^triai > £ler even outside this range, we infer that the 
LER energy is embedded inside a left half-circle in the 
complex energy plane: the left half of the circle where 
the exact solution is embedded according to Ref. [l8[ . We 
summarize it by saying that under the conditions stated 
above, the energy of the LER, Rei? LER = £ LER , which is 
associated with the real part of the complex eigenvalue 
of the non-hermitian hamiltonian, satisfies the following 
modified complex variational principle: 



minRe(*g|iJ e |* e ) 



(8) 



Hohenberg-Kohn theorem. Having established the 
plausibility of Eq.© for trial functions that can be ex- 
panded in terms of resonance wavefunctions, an analog 



But the opposite result is obtained by employing 
as a trial function to estimate the lowest-energy eigen- 
value of Hg^ . We conclude that the original assumption 

of Eq. © is impossible if v^ t and v^ t differ by more than 
a constant. 

To see the problem from a different perspective, we 
now examine the Levy-Lieb [Trl Il7j constrained search 
algorithm in the present context. The LER state is the 
one that, among all the normalized wavefunctions that 
make the complex energy (^fg\Hg\^fg) stationary, mini- 
mizes the expectation value of Re^glHg^g). Following 
Levy, we perform the minimization in two steps, first con- 
straining the search among all the wavefunctions yield- 
ing a prescribed complex density, {^g t— > rig}, and then 
among all possible complex densities, {ng}. The energy 
of the LER is then given by: 



minRe 



dz n e (z)v oxt (z) + Fg[n 9 ] 



s.t. {cl} , 
(10) 



Reify [rig] = min 



Re(^*g\Te + V ee ,g\Vg)\ s.t. {c2} 

(11) 



and constraints cl and c2 are as discussed before: 



cl : 



J dzn g (z) = N 



c2: —(9* e \H 9 \9 )=0 

OWg 



(12) 
(13) 



In spite of the formal resemblance of Eq. fTU]) with 
the density-variational principle that serves as a start- 
ing point to derive Kohn-Sham equations, condition c2 
makes of this a very different problem. It introduces a 
seemingly very complicated explicit dependence of Fg [ng] 
on v cx t 1 preventing proof of the HKT analog. 

But we now invoke Eq.©. According to it, constraint 
c2 can be lifted altogether. The resulting (unconstrained) 
search of Eq. fTTj) defines a universal functional Reify [ng], 
just as in the ground-state case. There is no explicit 
dependence of Reify [ng] on i> ox t, and the HKT analog is 
established. 
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To access the lifetime of the LER, denote by the 
wavefunction that, within the set of functions yielding ng, 
minimizes the real part of (*f?g\Hg\if?g) subject to the nor- 
malization constraint cl. It is a functional of ng, ^&g[ng]. 
If we further define IvaFg [ng] as the imaginary part of: 



Fg[ng] = {%[ng]\fg + V ee ,g \*g{ng]} 



(14) 



then the inverse lifetime of the LER is given by the imag- 
inary part of the sum of Fg[ng] and J G?zne(z)u e xt(z)- 

When apart from being the resonance of lowest energy, 
the LER is also the resonance of longest lifetime, a typ- 
ical case (e.g. our toy example), then Eq. (fT(j]) can be 
subsumed by a two-component minimization yielding at 
the same time the energy £ and inverse lifetime r = H/t 
of the LER: 



Re 
mm I OT 

ng \ — zlm 



J dz ng (z)ii cxt (z) + Fg[ng 

s.t.{cl} 



We have admittedly not addressed here the two fun- 
damental questions that immediately arise: (1) What is 
the best way to cast the complex analog of the Kohn- 
Sham scheme for practical calculations? and (2) What 
is the functional form of Fg[ng]l For one electron, it is 
simple to show that Fe[ng] = e~ F[ng], where F[ng] 
is the ground-state functional evaluated on the complex 
density. 

Our derivation applies to the LER of unbound systems 
such as negatively charged atoms or molecules. How- 
ever, using the Gel'fand-Levitan[zO] equation it is quite 
straightforward to extend our formulation to systems 
that support also bound states. It has been shown al- 
ready that using the Gel'fand Levitan equation one can 



remove bound states from the spectrum and obtain an 
effective potential which supports resonances onlv[2l|. 
However, from a numerical point of view it might be a 
heavy task problem since the computation of new effec- 
tive potentials that support the same resonances as the 
original problem, but not any of the N bound states, 
requires the often prohibitive calculation of those bound- 
state wavefunctions. Our extension of the HKT for 
the LER of unbound systems holds also for atoms and 
molecules in the presence of external DC or AC electric 
fields, since the field-free ground (bound) state becomes a 
resonance state as the DC or AC fields are turned on (for 
the calculation of such resonances via complex scaling see 
refs. and 0). 

Negative electron affinities. We comment briefly on the 
computation of negative electron affinities as measured 
1^ experimentally for many molecules via electron transmis- 
sion spectroscopy [22j]. The standard definition of the 
electron affinity is: A = E^ N ~^ - E^ where E^ N -^ is 
the ground-state energy of the neutral molecule, and E^ 
is the ground-state energy of the negative ion. The latter 
is precisely equal to E^- N ^^ when the ion is not bound, 
so A is zero in such cases. Confusion arises in practice 
when a finite basis set used in DFT calculations artifi- 
cially binds the ion and predicts finite (negative) values 
for A. But the experiments measure a different quantity: 
A = E^ N ^^ — Rel^ER and this is not directly accessible 
via standard DFT calculations in the limit of an infnitc 
basis-set. It is nonetheless interesting that A — A is ac- 
curately given in many instances by the error associated 
with the use of a finite basis set (see discussion in ref. [ 
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